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Abstract 

The gas pixel detector (GPD) dedicated for photoelectric X-ray polarimetry is selected as the focal plane detector for 
the ESA medium-class mission concept X-ray Imaging and Polarimetry Explorer (XIPE). Here we show the design, 
assembly, and preliminary test results of a small GPD for the purpose of gas mixture optimization needed for the phase 
A study of XIPE. The detector is assembled in house at Tsinghua University following a design by the INEN-Pisa 
group. The improved detector design results in a good uniformity for the electric field. Pilled with pure dimethyl 
ether (DME) at 0.8 atm, the measured energy resolution is 18% at 6 keV and inversely scales with the square root 
of the X-ray energy. The measured modulation factor is well consistent with that from simulation, up to ~0.6 above 
6 keV. The residual modulation is found to be 0.30% + 0.15% at 6 keV for the whole sensitive area, which can be 
translated into a systematic error of less than 1% for polarization measurement at a confidence level of 99%. The 
position resolution of the detector is about 80 jim in PWHM, consistent with previous studies and sufficient for XIPE 
requirements. 
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1. Introduction 

X-ray polarimetry is expected to be a powerful tool 
for astrophysics, offering extra information in addition 
to X-ray imaging, spectroscopy and timing. It allows us 
to probe the magnetic field via synchrotron radiation or 
test the geometry via scattering, and is capable of testing 
fundamental physics such as quantum electrodynamics 
and general relativity under extreme magnetism or grav¬ 
ity [1, 2, 3]. Despite the high demand in astrophysics. 
X-ray polarimetry has been an unexplored area for 40 
years since the experiments on the OSO-8 satellite [4] 
in 1970s, due to the absence of sensitive technology. 

Along with the development of micro-pattern gas de¬ 
tectors, it has become possible to image the tracks for 
electrons of a few keV in gas chambers, allowing for 
sensitive X-ray polarimetry depending on the photo¬ 
electric effect using the gas electron multiplier (GEM) 
with pixel readout [5, 6]. Compared with other readout 
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techniques, the gas pixel detector (GPD) offers symmet¬ 
ric measurement in the two dimensions and delivers a 
low systematic error below one per cent even without 
instrument spinning [6, 7, 8]. One of the key elements 
to the success of the GPD detector is the large-format 
high-resolution pixel readout chip. Several generations 
of dedicated ASIC chips have been developed, reach¬ 
ing to a pixel size of 50 jim and a chip size of 1.5cm 
X 1.5cm (105k pixels) [9, 10]. Two versions of sealed 
test chambers, a small one [11] and a large one [7], were 
designed by the INEN-Pisa group and assembled by Ox¬ 
ford Instrument Analytical Oy. The major improvement 
for the large version is that the background induced by 
the wall of the chamber is reduced and the electric field 
is more uniform [12]. 

In 2015, the X-ray Imaging and Polarimetry Explorer 
(XIPE) [3] was approved for phase A study by the Euro¬ 
pean Space Agency (ESA). It is a space telescope con¬ 
cept dedicated to X-ray polarimetry in response to the 
call for medium-class missions. One of our tasks is to 
build refillable sealed GPD detectors to test and opti¬ 
mize the gas mixture of the detector. We therefore de- 
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Figure 1: A schematic drawing of the GPD. 


scribe here the assembly of the GPD detectors and some 
preliminary test results. The optimization of the gas 
mixture and related tests of the detector will be reported 
in follow-up papers. Our work is based on the small 
version GPD, because it is relatively simple in structure 
and assembly, and sufficient for testing the gas mixture. 



Figure 2: Detector assembly, (a) ASIC chip in the ceramic package 
and both on the PCB; (b) GEM framing; (c) spacer and GEM; (d) 
gluing the spacer/GEM to the package; (e) titanium cap and beryllium 
window; (f) a completely assembled GPD. 


2. Detector structure and assembly 

We start with the design of the small GPD [6, 11], see 
Figure 1 for a schematic drawing. The ASIC chip [10] 
is mounted using low outgassing silver epoxy inside a 
vacuum tight 304-pin ceramic package, which is sol¬ 
dered onto a printed circuit board (PCB). The bottom 
of the chip is connected to the electronic background 
and the heat from the chip is dissipated into the base of 
the package and then to the PCB board, where a Peltier 
cooler can be mounted on the other side and multiple 
via holes are used to increase the thermal conductivity. 
The surface of the chip is measured to have an orien¬ 
tation tolerance of less than 20 jum with respect to the 
shoulder of the package using a measuring microscope. 
The distance between the chip surface and the package 
shoulder is 0.8 mm, which is the induction distance be¬ 
tween the anode and the bottom layer of the GEM. The 
GEM foil is manufactured by SciEnergy Inc. It is 50 
jum thick insulated by liquid crystal polymer; the laser- 
etched holes have a diameter of 30 pm and a pitch of 50 
pm in a hexagonal pattern. Eour knife-edge clamps are 
used to apply tensions to the GEM foil and a ER4 frame 
is fixed on it to keep the tension using double sided 
sticky tape. The framed GEM foil is then placed right 
above the shoulder of the ceramic package; the frame is 


just outside the package and is in the air side so that the 
choice of its material is not important. A ceramic spacer 
of 1 cm thick (the drift distance) stands above the GEM 
foil and supports the cap of the chamber, which is a ti¬ 
tanium plate with a 17mm x 17mm square hole sealed 
with a 100 pm thick beryllium window, both serving as 
the cathode of the detector. On the top of the titanium 
plate, a copper tube with a diameter of 6 mm is mounted 
for vacuum pumping and gas filling. All the parts for the 
main body of the chamber are glued together using low 
outgassing epoxy and cured at a temperature of 60 °C. 
The copper tube is welded onto the titanium plate and 
some epoxy is used to further seal the joint. After the 
chamber is filled with gas, the tube is cut and sealed us¬ 
ing ultrasonic welding. To change the gas mixture of the 
chamber, we will just cut the copper tube. Before filling 
the gas, the detector is pumped and baked for weeks un¬ 
til a vacuum down to ~10“® mbar. Some of the pictures 
during the assembly process and the final detector are 
shown in Eigure 2. 

Eor all the tests in this paper, we adopt a gas mix¬ 
ture of pure dimethyl ether (DME) at 0.8 atm, for a di¬ 
rect comparison with previous results [13]. The DME 
is purified to better than 0.99999 measured by gas chro¬ 
matography. 
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Figure 3: Inner structure of the chamber that connects the copper tube 
and the inner chamber, for the original design (left) and the improved 
design (right). 
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Figure 4: Intensity map of the electric potential in the middle plane of 
the chamber, for the original design (left; Figure 3a) and the improved 
design (right; Figure 3b). 

3. Uniformity of the electric field 

For the small GPD design, a major challenge is to 
improve the uniformity of the electric field, which will 
directly influence the electron drift and consequently the 
accuracy of the electron track measurement. In our first 
design (Figure 3a), a slot on the thick wall of the spacer 
and a right-angle hole on the titanium cap is machined 
to connect the copper tube and the inner chamber. How¬ 
ever, such a design creates an asymmetric structure and 
leads to large non-uniformity for the electric held. Fig¬ 
ure 4a shows the intensity map of the electric potential 
in the middle plane along the drift direction (Z - 5 mm 
above the GEM top plane) in the GPD chamber sim¬ 
ulated using Ansoft Maxwell, assuming a voltage of - 
2740 V on the drift plane and -740 V on the top plane 
of the GEM. As one can see, the electric held is quite 
uneven toward the tube side. To test the uniformity, we 
constructed a narrow beam of X-rays and illuminated 
the detector at different locations. At positions close to 
the tube side, the cumulative charge map produced by 
all the events is elongated, indicative of a distorted elec¬ 
tric held (Eigure 5a). The test did not cover the whole 
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Figure 5: Surface uniformity test using a narrow beam for the original 
design (top; Figure 3a) and the improved design (bottom; Figure 3b). 
The crosses indicate the incident positions for each beam. The con¬ 
tours are the collected charge distribution for each beam test, shown 
in log scale. 


sensitive area because the non-uniformity is already ob¬ 
vious and a full coverage takes a huge amount of time. 

To correct the held distortion, we changed the design 
of the interconnection structure. The spacer wall be¬ 
comes thinner and stands outside of the right-angle hole 
in the titanium plate (Eigure 3b). This makes the ce¬ 
ramic spacer look the same on each face from inside. 
The simulated potential map shows that the maximum 
differential voltage in the middle plane is only 2.4 volt 


3 









Table 1: The choice of Bragg crystals and diffraction orders, corre¬ 
sponding X-ray energies, and measured modulation factors. 


E (keV) 

crystal 

order 

F(%) 

2.67 

MgF2 

1st 

27.2+1.3 

3.74 

Al 

1st 

42.3+0.4 

5.33 

MgF2 

2nd 

54.1+0.6 

6.09 

LiF 

2nd 

58.2+0.7 

7.49 

Al 

2nd 

61.9+1.3 


Table 2: Spectral fitting results at different energies. 


E 

(keV) 

Peak 

(ADC) 

FWHM 

(ADC) 

AEIE 

(%) 

2.67 

1901.1 + 2.3 

230.0+ 1.8 

26.10 + 0.23 

3.74 

2744.9 + 0.9 

278.2 + 0.7 

22.44 + 0.06 

5.33 

3938.7+ 1.8 

331.7+ 1.4 

18.97 + 0.08 

6.09 

4548.6 + 2.3 

358.1 + 1.9 

17.83 + 0.10 

7.49 

5650.3 + 5.4 

402.0 + 4.9 

16.23 + 0.21 


or 0.14% (Figure 4b). The narrow beam test indicates 
that the charge distribution is no longer elongated or po¬ 
sition sensitive (Figure 5 a). The uniformity of the elec¬ 
tric field is much improved and such a design is capable 
of testing the performance of the detector over the whole 
sensitive area. 

4. Test and results 
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Figure 7: Effective gain of the GEM versus the GEM voltage. 



100 200 300 400 500 600 

Effective Gain 


In order to characterize the detector performance, we 
constructed a desktop X-ray source including an X-ray 
tube manufactured by Oxford Instrument and crystals 
to produce fully polarized monochromatic X-ray beams 
by Bragg diffraction at 45 degrees, following the setup 
described by Muleri et al. (2008,2010) [14, 13]. The 
choice of crystals and their Bragg diffracted energies at 
45 degrees are listed in Table 1. The bremsstrahlung 
continuum from the X-ray tube serves as the incident 
beam for diffraction. A capillary plate manufactured 
by Flamamatsu, which has a hole diameter to thick¬ 
ness ratio of 1/167, is mounted in the light path after 
the diffraction to select photons due to diffraction rather 
than Compton scattering. With such a source, we are 
able to produce polarized X-rays at 2.67, 3.74, 5.33, 
6.09, and 7.49 keV. 

4.1. Spectral response 

The incident X-rays have a narrow beam size due 
to collimation with the capillary plate. The beam is 


Figure 8: Fractional spectral resolution (FWHM/E) at 3.74 keV ver¬ 
sus the effective gain of the GEM. 

aligned to illuminate the central region of the detector. 
The analysis of the data is implemented in the PIXY 
package developed by the INFN-Pisa group. We select 
events in a circular aperture with a radius of 1 mm. A 
noise threshold of 8 ADC is adopted. The energy for 
each event is calculated by summing all of the pixel val¬ 
ues after the noise cut. The energy spectra are shown in 
Figure 6. The photopeaks due to the 1st and 2nd orders 
of diffraction with the MgF 2 and Al crystals are well 
distinguished. The energy spectra also show a tiny low 
energy tail (e.g., see the 6.09 keV spectrum), which is 
more prominent for high energy photons and is caused 
by absorption of the photoelectron on the GEM foil or 
the beryllium window that terminates further ionization 
in the gas [7]. To extract the spectral information, we fit 
a Gaussian function to each peak with results listed in 
Table 2. 


4 










3500 



8000 

3.74 keV 

2000 

6.09 keV 

3000 


5.33 keV 


I 



2500 

- 

i 

6000 

1 

1500 

1 

2000 

- 

1 

a 

1 


1 



■ 

1 4000 

M 

1 1000 


1500 


■ 

u 

■ 



1000 

i j 

1 

2000 

1 

500 

A 

500 

k i 

m 


7.49 keV 



0 

^ J 


0 


0 

_ 


0 2000 4000 6000 8000 0 2000 4000 6000 8000 0 2000 4000 6000 8000 

Pulse Height (ADC) Pulse Height (ADC) Pulse Height (ADC) 


Figure 6: Pulse height distribution measured with X-rays created by 45-degree Bragg diffraction with MgF 2 (left), A1 (middle) and LiF (right) 
crystals. 



Figure 9: The pulse hight of the photopeak versus the incident X-ray 
energy. 

The drift field is fixed at 2 kV/cm, where the coeffi¬ 
cient of transversal diffusion for electron drifting is min¬ 
imized in pure DME at 0.8 atm. The induction field (or 
the collection field) is set to 3 kV/cm; a smaller field 
may lower the fraction of electrons collected by the an¬ 
ode, and a higher field may lead to multiplication in the 
induction region. Figure 7 shows the change of the ef¬ 
fective gain as a function of the GEM voltage, and Fig¬ 
ure 8 illustrates the variation of the spectral resolution 
along with the change of the effective gain. We choose 
a GEM voltage of 500 V where the effective gain is near 
400, and the spectral resolution is minimum and insen¬ 
sitive to the change of the gain. We note that the choice 
of this GEM voltage is also optimal for polarimetry, re¬ 
sulting in a maximum degree of modulation, which is 
not shown here. 

The detector exhibits a good linearity between the in¬ 
cident X-ray energy and the output pulse height (see 
Figure 9). The measured energy resolution (Figure 10) 
in terms of FWHM/E is found to be scaled with the 
square root of the X-ray energy, consistent with that the 
noise is dominated by Fano fluctuation. 



Energy (keV) 

Figure 10: Measured energy resolution (AE = FWHM) versus the X- 
ray energy, along with a best-fit curve in the form of AE/E oc 1/ ^fE. 


4.2. Polarimetry 

Examples of electron tracks produced by X-rays of 
different energies are displayed in Figure 11. For po- 
larimetric analysis, events with two or more separate 
clusters of charges are not used. Those having a total 
hit of less than 27 pixels are removed too, as they may 
cause a systematic effect at +60 degrees due to the pat¬ 
tern of the GEM holes. For a direct comparison with 
previous results, we further discard 25% of events with 
low eccentricity following Muleri et al. (2010) [13]. The 
emission angle of photoelectrons is derived by finding 
the principal axis of the whole track for events with en¬ 
ergies below 3 keV, or using the so-called impact point 
method [15] for high energy events. 

A modulation curve, which is the histogram of pho¬ 
toelectron emission angles, measured with a fully po¬ 
larized beam of 3.74 keV is shown in Figure 12 . 
The degree of modulation is defined as the fractional 
peak-peak amplitude of the modulation curve, (Nmax - 
iVmin)/(Vmax + Nmin)- The modulation factor jj. is the 
degree of modulation in response to a fully polarized 
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Figure 11: Example electron tracks produced by X-rays of dilferent 
energies. 


beam and is one of the most important parameters for a 
polarimeter, which is directly scaled with the sensitiv¬ 
ity for polarimetry. The measured modulation factors 
at all of the available energies are listed in Table 1 and 
displayed in Figure 13. For comparison, the results ob¬ 
tained by Muleri et al. (2010) with the same gas mix¬ 
ture are also shown [13]. The curve in the figure repre¬ 
sents the expected modulation factor versus energy cal¬ 
culated using GEANT4/Garfield simulations taking into 
account the transversal electron diffusion during drift, 
multiplication and collection, electronics noise, and the 
resolution of the readout chip. Our result is well consis¬ 
tent with that obtained with the simulation. 

To test the systematics of the detector, we constructed 
an unpolarized source by illuminating a manganese tar¬ 
get using the same X-ray tube, which produces fluores¬ 
cent lines at 5.9 keV. To further smooth out any possi¬ 
ble polarized component in the X-ray beam, e.g., due 
to scattering, we kept rotating the source slowly around 
its beam axis during the test by 1 degree per second for 
some 80 complete cycles. About 960,000 valid events 



<p (rad) 

Figure 12: Modulation curve measured with a fully polarized beam of 
3.74 keV. The solid line is the best-fit cosine function with a degree of 
modulation of 42.3% ± 0.4%. 
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Figure 13: Modulation factor versus energy. The red circles indicate 
results measured with our detector, while the blue diamonds are re¬ 
sults obtained by a previous detector with the same gas mixture [13]. 
The solid curve is the result obtained using GEANT4/Garfield simu¬ 
lations. 


were detected by the detector, almost evenly distributed 
on the sensitive area, 15mm x 15mm. This, given a 
modulation factor of 0.58 at 5.9 keV, corresponds to a 
minimum detectable polarization (MDP) of 0.75% at 
99% confidence level. Figure 14 shows the measured 
modulation curve for the unpolarized beam. The mea¬ 
sured degree of modulation is 0.30% + 0.15%, corre¬ 
sponding to a 99% upper limit of 0.55% in degree of 
modulation or 0.94% in degree of polarization by taking 
into account the modulation factor. The result is consis¬ 
tent with statistical fluctuation from zero polarization, 
and indicates that the instrument systematic error is be¬ 
low 1% for polarimetric measurement. 

4.3. Imaging 

We constructed a narrow X-ray beam to characterize 
the imaging properties of the GPD. The X-ray tube has 
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Figure 14: Residual modulation measured with an unpolarized beam 
at 5.9 keV. The degree of modulation is 0.30% ±0.15%, coiTesponding 
to a 99% upper limit of 0.55% in degree of modulation or 0.94% in 
degree of polarization. The bottom panel is a zoom-in of the top panel. 


a focal spot size of about 100 yum. We use a stainless 
steel diaphragm with a diameter of 30 /rm mounted at 
the end of a 18 cm long tube to collimate the source. 
The beam size is measured by scanning the beam with a 
Si-Pin detector at the same distance of the GPD, which 
is half masked with a sharp edge [16]. Assuming that 
the intensity distribution of the X-ray beam can be de¬ 
scribed by a 2D Gaussian, an error function is adopted 
to fit the variation of the counting rate versus the scan¬ 
ning position. The beam size in FWHM is measured to 
be 31.8 jum in one direction and 28.6 pm in the orthogo¬ 
nal direction, consistent with the size of the diaphragm. 

Figure 15 shows the detected image by the GPD with 
the narrow beam. The X-ray interaction point is adopted 
for image reconstruction. A 2D Gaussian function is ad¬ 
equate to fit the data with;t'^ = 396.7 and 393 degrees 
of freedom. The measured FWHM of the point spread 
function is 87.4+0.9 pm in the X direction and 79.8+0.9 
pm in the Y direction. Taking into account the contri¬ 
bution of the beam size, the intrinsic spatial resolution 
is about 78 pm in FWHM. We take three exposures by 
shifting the beam by 300 pm in two orthogonal direc¬ 
tions. Figure 16 shows an image by combining all of 
the data, to simulate an image of three nearby point-like 
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Figure 15: Image of a narrow X-ray beam detected by the GPD and 
fitting with a 2D Gaussian profile. The beam size is around 30 ^m in 
FWHM and the intrinsic point spread function is measured to have a 
FWHM of about 78 fim. 



Figure 16: Image of three sources with 300 fim apart in X and Y. 


sources. 

5. Discussion and Conclusion 

In this paper, we presented the assembly of the GPD 
and some preliminary test results. We demonstrate that 
we are able to assemble detectors qualified for the phase 
A study of the XIPE mission. The improved design of 
the detector structure results in a good uniformity for the 
electric field. Both the simulation with Maxwell and test 
using a narrow X-ray beam suggest that the electric field 
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for the improved design is capable of X-ray polarimetry 
with low systematics. To compare with previous results, 
mainly with Muleri et al. (2010) [13], we choose a gas 
mixture of 0.8 atm pure DME that is relatively easy to 
fill (no need to mix with other gas components). The 
GEM foils used in the two studies are different. Our 
GEM is 50 yum thick with holes of 30 /rm in diameter 
and 50 /tm pitch, while their GEM is 100 jum thick with 
holes of 48 //m in diameter and 80 /rm pitch. The choice 
of drift field, GEM voltage, and the induction field is 
to maximize the imaging precision and signal to noise 
ratios. 

The spectral resolution of our detector at 6 keV is 
better than that reported in Muleri et al. (2010), 18% 
versus 24% [13], which may be a direct consequence of 
the impurity of the DME gas. We once filled a detector 
and obtained worse resolution, but it was improved after 
we resolved the outgassing problem of the gas filling 
tubes between the gas cylinder and the detector. They 
also obtained nice spectral resolutions, e.g., 15% [11] or 
19% [7] at 6 keV with a different gas mixture (DME:He 
8:2 at 1 bar). Thus, the previous low resolution with 
pure DME is not due to the nature of the gas, rather it 
could be due to the non-uniformity of the GEM gain 
at different positions, as their beam size is much larger 
than ours. Once filled and sealed, our detector does not 
show a decline of performance with time, at least on a 
timescale of a year or so, which will be reported in detail 
later. 

The modulation factors that we measured is slightly 
higher than in Muleri et al. (2010) (Eigure 13). We note 
that the impurity of gas, which leads to electron attach¬ 
ment and loss during drifting, is not crucial to the mod¬ 
ulation measurement. We found that consistent modu¬ 
lation factors could be obtained with the same detector 
filled with gases of different purity, e.g., with an energy 
resolution of 18% to ~30%. We note that the difference 
on the GEM thickness and hole resolution by the two 
studies, according to simulation, has almost no effect 
on the modulation factor with the choice of such a gas 
mixture. 

The residual modulation, measured to be 0.30% + 
0.15%, is consistent with the value 0.54% obtained by 
Muleri et al. (2012) [7] and 0.18% + 0.14% as reported 
by Bellazzini et al. (2013) [8], although the two previ¬ 
ous studies used a different gas mixture (DME:He 8:2 
at 1 bar). The measurement by Muleri et al. (2012) 
[7] should be dominated by statistical errors, because 
the number of counts we collected is much higher than 
what they got (960k versus 125k). Limited by the source 
intensity, our measurement for the residual modulation 
covers the whole sensitive area of the detector. Oth¬ 


erwise, it would take too much time to get a mean¬ 
ingful measurement just in the focal point region us¬ 
ing our existing setup. If we extract events in the cen¬ 
tral region of 3mm in radius, the residual modulation is 
0.15% + 0.39% constructed from 129k counts, which is 
dominated by statistical errors. 

The imaging performance of our detector is well con¬ 
sistent with that measured by Soffitta et al. (2013) [16]. 
A position resolution of about 80 //m in EWHM cor¬ 
responds to an angular resolution of 4.7 arcsec given a 
telescope with a focal length of 3.5 m, sufficient to sam¬ 
ple the point spread function for the XIPE optics which 
is required to have a angular resolution of less than 30 
arcsec. 

To conclude, our assembled GPD shows a perfor¬ 
mance consistent with that of the best ones in previous 
studies, which will enable us to test and optimize the 
gas mixture for the phase A study of XIPE. 
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